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Recent results for the finite-temperature static and dy- 
namical properties of the planar t — J model, obtained by 
a novel numerical method for small correlated systems, are 
reviewed. Of particular interest are T > charge and spin 
response functions: optical conductivity and spin sus- 

ceptibility x(9, i^)) which show very universal features at in- 
termediate doping. In spin dynamics the universality is best 
seen in the local spin correlation function, which appears to 
be T- independent and only weakly doping dependent. Re- 
sults apply directly to the inelastic neutron scattering and to 
the NMR relaxation experiments in cuprates. Analysing (t{u)) 
we find that the current correlation function C{ijj) is nearly 
T- and moreover tj-independent in the regime T < J, leading 
to the (j{ui) behavior, very close to the marginal Fermi liquid 
concept and experiments. These properties are interpreted in 
connection with the large entropy, related to the large degen- 
eracy of the low-lying states and quantum frustration in the 
doped antiferromagnets. On the other hand, results for the 
uniform susceptibility Xo{T), as well as for the cr{oj), reveal 
qualitative changes on entering the underdoped regime. 



I. INTRODUCTION 

Strongly correlated electrons have been in recent years 
intensively investigated theoretically, mainly in connec- 
tion with the open question of the mechanism of super- 
conductivity as well as the unusual metal properties in 
cuprates. Here we will discuss only the anomalous nor- 
mal state properties of cuprates at T > T^. These are 
generally attributed to correlation effects, nevertheless a 
unifying theory still remains a challenge. 

In this article we will address several questions re- 
lated to the theoretical understanding of the normal-state 
properties of cuprates: 

a) Can strong correlations alone account for the unusual 
normal-state properties, as manifested by the charge re- 
sponse: optical conductivity cr(a;) and d.c. resistivity 
Po, and the spin response: dynamical spin susceptibility 
x{q,<jj), related to the NMR relaxation Ti, etc.? 

b) What is the origin of the unconventional universality 
within the 'normal' metallic state of cuprates, in partic- 
ular in the 'optimum-doping' regime, obtained by doping 
the reference antiferromagnetic (AFM) insulator with the 
mobile hole-type charge carriers? 

c) How many relevant energy scales exist in the doped 
AFM? 

d) Are the simplest prototype models, such as the t — J 
model (and analogous Hubbard model), enough to re- 
produce qualitatively and quantitatively experimental re- 
sults in cuprates? 



f) Can we learn something by studying small correlated 
systems at T > 0? 

From the recent experiments testing the normal-state 
properties such as the d.c. resistivity Q], optical conduc- 
tivity HJ^, neutron scattering UJ|], and the NMR and 
NQR relaxation |^], a unifying picture seems to emerge 
[Q. These properties mainly depend on the extent of 
doping, and the regimes have been classified into the un- 
derdoped, the 'optimal' (intermediate-doping), and the 
overdoped, respectively. Whereas the 'optimal' doping 
refers in a narrow sense usually to samples with the high- 
est Tc, the same broader region of doping, with effective 
hole concentration approx. Ch = 0.15 — 0.25 also 
exhibits most universal properties. Down to the lowest 
T resistivity follows the remarkable po oc T law, the 
uniform susceptibility xo is varying monotonously with 
T, etc. On the other hand, the underdoped materials (ef- 
fective Ch ^ 0.15) are characterized by the onset of the 
'pseudogap' at T < T 3> T^, which shows up in the kinks 
in p(T), xoiT) etc. 

We investigate in this review mainly the intermediate- 
doping (optimal) regime. We discuss finite-T (normal- 
state) properties of the doped AFM within the t — J 
model 
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where c|^(cis) are projected fermionic operators, pro- 
hibiting double occupancy of sites. The ground state 
of the t — J model and related properties have been in- 
tensively studied both by analytical Q] and numerical 
methods nevertheless some basic questions, e.g. con- 
cerning possible pairing in the ground state and the origin 
of anomalous metallic properties at finite T > 0, remain 
unsolved. 



II. FERMI LIQUIDS 

Normal metals are expected to follow the concept of 
the Landau Fermi liquid. The excitations in this case 
should be quasiparticles with well defined q and energy 

near the Fermi energy, with the damping T{q)~^ oc 
(fig- — EfY- Consequences are well known, e.g. for the 



dynamical susceptibility 
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leading to the Korringa law NMR ^ cx T. Also, the 
Drude form is expected for a{u!) with cx etc. 

The normal-state properties of cuprates, however, do 
not follow the Landau scenario. Experiments have been 
mainly described within two alternative concepts. The 
nearly antiferromagnetic Fermi liquid has been in- 
troduced by Moriya and Millis et al. |^ to describe the 
spin fluctuations in the doped AFM, mainly in connec- 
tion with the anomalous NMR relaxation phenomena 
in cuprates. The low-w spin response is modeled with 
the overdamped spin fluctuations peaked at the AFM 
wavevector Q = (tTjTt), 



(3) 



AFM correlation length ^(T) is assumed to be critical 
in the doped regime, in particular with C^(r) cx 1/T 
in the z = 2 dynamical RPA-type theory while the 
postulated mapping to the quantum critical regime of the 
nonlinear sigma model would imply ^(T) (x 1/T |l0| . 

Alternative scenario has been summarized within the 
marginal Fermi liquid (MFL) hypothesis Here, 
the spin and charge susceptibility have been postulated 
in the anomalous form 



x"{q,Lo) ^ C Sgnto, \u;\>T, 



(4) 



leading to the marginality of the quasiparticle self energy 

ttT — iuj 



J:{lu) - X[2uj In- 
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Consistent with Eq.(||), the low-frequency cr(w) has been 
written in the generalized Drude form introducing an 
anomalous effective relaxation rate t^^ ^ 2ttX{uj + ttT). 



III. FINITE TEMPERATURE METHOD 

r > properties of the t — J model have been so far 
calculated only by the high-temperature series expansion 
and via full or partial diagonalization of quite small 
systems. Recently the present authors introduced a novel 
numerical method, based on the Lanczos diagonalization 
method combined with random sampling [p^ . The lat- 
ter method has the advantage that it allows the study of 
dynamical and static response functions within the most 
challenging regime T,lo < J, where one could hope for 
a universal behavior, as observed in the experiments on 
cuprates. It should be noted that analogous Hubbard 
model seems more feasible for T > studies, due to the 
possible application of the Quantum Monte Carlo meth- 
ods 1^. Nevertheless, in the doped case the well known 
minus sign problem so far mainly prevented calculations 
within the regime T,uj < J H . 



Let us present briefly our method for T > calcula- 
tions, showing in more detail the evaluation of the statical 
expectation value of an operator A, 
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where /? — 1/T (we use = h = 1 furtheron) and the 
sum runs over the chosen complete basis set of orthonor- 
mal wavefunctions \n),n — 1,N, spanning the Hamilto- 
nian H. If we could perform the full diagonalization of 
the problem and find all eigenstates l^";) and correspond- 
ing energies we would express the result in a usual 
way, 
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We choose an alternative approach. With each basis 
function \n) we start a Lanczos procedure \4)q) = 
generating an orthonormal set of functions |(/)J^),to = 
0,M, 
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For the chosen number of Lanczos steps M <^ N we 
then diagonalize the tridiagonal matrix of coefficients 
bnm to find the energies and the corresponding 



eigenfunctions iV'm)- Within the restricted basis of these 
functions one can write the approximation for Eq.(^ as 

No M 



No M 
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It is evident that the average evaluated via Eq.(||) is 
equivalent to Eq.(0) for the sampling over the full ba- 
sis set Nq^ N and M ^ N - 1. 

Our claim is that very accurate results can be obtained 
via Eq. even for a severly reduced number of Lanczos 
steps M <^ N and for a partial random sampling of the 
basis states Nq ^ N, instead of the full sampling. Let us 
first investigate the method for the full sampling Nq — N. 
We note that the Lanczos procedure Eq. represents an 
iterative action of the operator H on the initial function 
\n). Performing the (high temperature) expansion in /? 
of the numerator and of the denominator of Eq. (j^) , it is 
easy to prove that both power series are correct up to the 
order M. On the other hand, (A) evaluated via Eq.(|^) 
is accurate also for T = (provided M > Mq), since 
it is a standard experience that the Lanczos procedure 
converges quite rapidly, typically in M = Mq ~ 50 <C 
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steps, from an arbitrary initial function |n) (not orthog- 
onal to the ground state) to the ground state energy Eq 
and the corresponding wavefunction |^'o)- Hence, at full 
sampling our Lanczos-based method represents an opti- 
mum interpolation scheme between the high temperature 
T ^ oo and the ground state T — result. 

The second important ingredient is the reduction of the 
sampling to the partial random one with Nq <^ N. This 
step cannot be justified rigorously, but seems plausible in 
analogy with the statistical Monte Carlo methods. The 
most severe test is expected to be the result at T = 0, 
which can be anyhow tested via an alternative direct eval- 
uation with the known |^'o). The general experience is 
that good convergence is achieved for Nq <^ N. 

Even more appealing is the application of the method 
to the calculation of dynamical quantities. We choose 
as an example the optical conductivity, expressed within 
the linear response theory as 

cr{cj) = Re / dt e''^'C(t), 

^ Jo 

C{t) = mj) = 5](n|e(-''+^*)^je-^^*j-|n). (10) 

n 

In order to get the result at arbitrary /3 a double sum 
over all eigenstates is required. Instead, we 

use the approach described above for {A) generalized to 
dynamical C(t), 

No M 

C{t) (n|C>e-^^'""e^(^""-'~-^* x 

n m,k 

mMk)mj\n)- (11) 

Here, \tpk) f-i^d ^nk are generated by the Lanczos proce- 
dure analogous to Eq.(g) via the orthonormal basis set 
\(t>^), but with the initial condition 

l^o>-jl">/yWR- (12) 

While for Nq — M = N the expression (jl^) is equivalent 
to the exact one, the arguments for using Nq <C N and 
M <^ N are analogous to those described above for the 
static case. E.g., for the full sampling Nq ~ N Eq.(pT|) 
yields correct series in — /3 -f it and it, respectively, up to 
the M'^ order. Accordingly, (j{u!) has correct frequency 
moments (ujP), p — 0,M, for /3 ^ etc. Again for 
M > Mq the method yields meaningful result also for 
T — 0, which can be tested via the usual Lanczos method 
for the T — dynamical quantities § . 

It is important to stress that computational require- 
ments of the method for a system of the given size are 
comparable to those for the ground state evalutions via 
the Lanczos method, whereby the CPU time clearly in- 
creases with the number of random samples Nq. The 
CPU time and memory requirements are also enhanced 



due to the necessary reorthogonalization of Lanczos func- 
tions and due to the evaluation of matrix elements in 
Eq.(|ll). 

Let us finally comment on the finite size effects, which 
due to the restricted system sizes still remain the main 
obstacle to achieve reliable results for low T. It is a plau- 
sible observation that the method yields macroscopic- 
like results for T > T*, where T* is roughly determined 
by the average level separation in the low-energy many- 
body spectrum, being dependent on the system size. For 
T < T* the finite-size effects become pronounced, e.g. 
dynamical quantities in general reveal (spurious) peak 
structures etc. It is however very important to note that 
within the t — J model, the T* is smallest within the in- 
termediate (optimum) doping regime 2/16 < Ch < 4/16, 
where we typically reach T* ^ 0.1 t for systems N = 
16, 18. The T* is larger both within the underdoped and 
the overdoped region (for the fixed system size). 

IV. SPIN SUSCEPTIBILITY 

Let us consider the dynamical spin response, as given 
by the susceptibility xiQ^ '-^) ^^'^ the corresponding dy- 
namical spin correlation function S{q, u) 

x"(<f,c.)-(l-e-^-)S(q,c.), 

S{q, cu) = Re dt e'"^* (S^(t)Slq~) . (13) 
Jo 

It has been shown by the present authors ||l^ that 
x"{q,Lo) shows the coexistence of the high- frequency 
[lj cx t) free-fermion-like contribution and the low-w spin- 
fluctuation contribution. At the same time, x"((f, to < T) 
reveals a pronounced T dependence, consistent with the 
MFL form, and (even quantitatively) with the NMR re- 
laxation in cuprates. 

It appears quite helpful to concentrate on the local 
spin correlation function S'l(w) p6[ | and its symmetric 
part S{uj) 

S{uj) = Sl{uj) + Sl{-uj) = (1 + e-P^)SL{uj). (14) 

It should be pointed out that Sl{^) and the related sus- 
ceptibility xl(^^) are directly measured (in cuprates) by 
neutron scattering [^||. The NMR relaxation as well 
yields the information on Sl{uj 0) provided that 
the AFM spin fluctuations q^Q— (tt, tt) are dominant. 
An important restriction for S{ll!) is the sum rule 

poo _ 

S{u;)du; = 7T{{S!f)^-{l-c,), (15) 
where Ch = N^/N is the hole concentration. 
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We perform the evaluation of S(lu) via Eq.(|lJ) by 
calculating S{q,u!), using the finite-T diagonalization 
method for small systems, in this case for the t — J 
model on the square lattice with N = 16, 20 sites. We 
fix J/t = 0.3 to remain in the regime of cuprates 0. 
We stress again that the method yields macroscopic-like 
results for T > T* > 0.1 t. _ 

In Fig. 1 we display the S{iLj) for Ch = 1/20,3/16 and 
several T in the range 0.1 < T/t < 0.7 It is imme- 
diately evident that S{uj) at 'optimal' doping Ch — 3/16 
is essentially T-independent in a wide T-range, although 
one crosses the exchange-energy scale T ^ J. For the 
underdoped case Ch = 1/20 the behavior is analogous 
for higher T > Tq ^ 0.7 J, consistent with the quan- 
tum critical regime within the AFM [0. Deviations at 
lower T < Tq (where the renormalized classical regime 
[0 is expected in the AFM) could be an indication for 
the onset of a 'pseudogap'. 



1.2 




FIG. 1. Local spin correlation function S{iij) for Ch=l/20, 
3/16 and various T. 

To follow the doping dependence we present in Fig. 2 
the variation with Ch at fixed T = 0.2 t < J, plotting the 
integrated intensity 



S{uj')duj' . 



(16) 



Again, we notice that for the chosen T the results are 
most reliable at the intermediate doping. The most 
striking message is that the initial slope of /s(w) and 
consequently S'l(ijJ ~+ 0) are nearly doping indepen- 
dent for < c/i < 0.25. This is well consistent with 
the NMR (NQR) relaxation rates T^^ as measured in 
La2-xSrxCu04 in the range x — ~ 0.15 |||,|l^. 

Only for the overdoped systems with Ch > 0.25 the 
low-frequency behavior changes qualitatively, where the 
latter part is strongly suppressed as expected in (more) 



normal Fermi liquids. /s(a; > J) is doping dependent 
even for < 0.25, consistent with the c^-dependence 
of the sum rule, Eq. (|l5|) . In addition, at the intermedi- 
ate doping S{ui) decreases smoothly (see Fig. 1) up to 
uj ~ At, this being the consequence of a free-fermion-like 
component |jl^. On the other hand, in the underdoped 
regime the dynamics is restricted to w < 3 J < t. 
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FIG. 2. Integrated spectra Is (a;) at fixed T = 0.2 t and 
various Ch- 

How can one explain the universality of S{u!) at in- 
termediate c/i? First, we note that up to Ch ^ 0.3 the 
dominant scale of spin fluctuations remains related to J. 
From the explicit expression in terms of the eigenstates 
of the system 



S{u;) = (1-f e-'=*-)|^e-'3^"|(n|Sr|m)p x 

n.m 

S{lj~E^ + E^), (17) 

one would conjecture (quite generally) a T-independence 
of response for lo ^ T. Although such plausibility ar- 
guments have been used previously, their validity clearly 
depends on the character and on the density of low-lying 
many-body states. To explain also the T-independence 
S{ijj < T), we only need to recall the sum rule 
Eq. (|l5|) and to assume that there is no characteristic scale 
LUc < T, which could introduce an additional low-w struc- 
ture in S{lu). 

A natural scale for the AFM is the (gap) frequency 
~ c/^, where c is the spin wave velocity and ^ the 
AFM correlation length. Here originates the essential dif- 
ference between the undoped and the 'optimally' doped 
AFM. While for an AFM in the renormalized classical 
regime ^ is exponentially large for T <C J, and con- 
sequently ujc < T [0, in the doped case ^ < l/y^ 
is determined predominantly by c/j, so ^ is rather T- 
independent for T < J, excluding oj^ < T < J. In other 
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words, the observed anomalous low-w behavior within the 
t — J model does not originate in a strong T-dependence 
of ^, as conjectured by the nearly AFM Fermi liquid, 
Eq.(|. 

We conclude the discussion of spin dynamics by the 
consequences of the universality of S{ui). The local sus- 
ceptibility is given by 

x'/M=tanh(^)sH. (18) 

Since neutron scattering probes only uj < J , one can 
simplify Eq.(p^) further by S{ljj) ^ Sq. Such form has 
been recently used to describe experiments It is also 
qualitatively consistent (but as well more restricted in 
the form) with the MFL Ansatz, Eq.(|) Q. 
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V. OPTICAL CONDUCTIVITY 

Let us investigate in an analogous way the dynamical 
conductivity cr(w), Eq. (|lo|) . Results for this quantity by 
the present authors [Q, obtained by the finite-T method, 
can be summarized as follows: a) in the intermediate- 
doping regime <j(lS) shows a non-Drude fall-off, consis- 
tent with the (jj-dependent relaxation rate within the 
MFL concept Eq.(||) deduced from experiments in 
cuprates and b) qualitative as well as quantitative 

results for cr(w) and p(r) agree well with the experimen- 
tal ones. 

Unlike S'(a;), C(w) does not obey a T-independent sum 
rule. Nevertheless, motivated by the universal spin dy- 
namics we reexamine in an analogous manner our results 
on (T(a;) [Q. In Fig. 3 we present the corresponding in- 
tegrated spectra /c(w) = /q" C{io')duj' at fixed doping in 
the 'optimal' regime for various T < t @. We establish 
that for T < J spectra Ic{^) are essentially indepen- 
dent of T, at least for the available T > T* . At the 
same time the slope of Ici^ < 2 t) is nearly constant, 
i.e. C{uj) ~ Co in a wide w-range, Cq being weakly J- 
dependent (tested for J/t = 0.2,0.6). 

We note that such C{uj) ^ Cq implies a nonanalytic 
behavior of a{Lu — s- 0), starting with a finite slope at 
w = 0. Moreover, with C{uj) = Co we can claim a simple 
universal form for co < 2 t 

a{Lo) = Co . (19) 

At first sight, it appears rather surprising that this form 
can be well fitted (for uj,T <C t) with a Drude-type 
form with an MFL- like effective relaxation rate = 
2nX{uj + riT) with specific A ~ 0.09 and ry - 2.7 [|l|]. The 
form Eq.(|lS|) for a{uj) trivially reproduces the remarkable 
linear law p oc T in cuprates ||l| , as well as the non-Drude 
falloff at w > T. 



FIG. 3. Integrated current-correlation spectra Iciy) for 
Ch =3/16 and various T. 

A quantitative comparison (assuming J/t — 0.3 and 
t — 0.4 eV) of <j{uj) with the experimental results 
is quite remarkable |]l8|, realizing the simplicity of the 
model. It is however evident that the expression ( p^ ) is 
more restricted than the MFL Ansatz, hence it should be 
retested in more detail with experiments. 

It should be mentioned that C(w) ^ Cq has been de- 
rived for a single hole conductivity within the retraceable 
path approximation jl9|, with a restricted validity for 
T > t (or possibly T > J). We find this behavior only 
for the intermediate doping, hence new arguments are 
needed. We can follow the analysis analogous to S{uj), 
Eq. (|l7|) , expressing C{oj) in terms of the eigenstates. As 
before, the T-independence of C{u! > T) seems plausible. 
The fact that C(w < 2t) ~ Co requires, however, that the 
current relaxation is very fast, i.e. determined only by the 
incoherent hopping and by the interhole collisions. The 
spin fluctuation scale J does not enter directly, i.e. even 
for T <^ J the spin system only serves as a random bath 
for the charge degrees of freedom (holons). This con- 
clusion remains valid as far as there is no characteristic 
frequency lOc <T (e.g. the 'pseudogap') in the system. 

On the other hand, our results also indicate the pos- 
sible qualitative changes when entering the underdoped 
regime. E.g. the cu = 2/18 doping seems to repre- 
sent a crossover to the latter behavior, where in contrast 
(e.g. for Ch = 1/20) we find a pronounced T- and lu- 
dependence of C{uj) at T < J. This qualitative change 
can be again attributed to a 'pseudogap' appearing at 
larger T in the underdoped systems ||l| . 

VI. ENTROPY 

Above arguments, both for Sl{uj) and G{uj), require 
a large density (degeneracy) of the low lying many-body 
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states, apparently being a crucial feature for the most 
challenging intermediate-doping regime. To quantify this 
statement we calculate the entropy density s = S/N, 
using again the finite-T diagonalization method, being 
less space and time consuming for static quantities [ p^ . 
In Fig. 4 we present the results for s{ch), obtained for 
iV 16 and iV = 20 at different T. The main lesson 
from Fig. 4 is that the 'optimal' cases with Ch ^ 0.15 — 0.3 
are characterized by the largest entropy s at low T < J, 
e.g. s > 0.2/site for T ~ 0.2 J, being almost one half 
of s{T = oo) for an AFM. This implies a very large 
degeneracy of the low-lying states, which could be at- 
tributed to the spin subsystem, frustrated by the hole 
motion. Calculated entropy, large even at moderate T 
and strongly dependent on doping, seems to be (also 
quantitatively) consistent with the specific heat measure- 
ments on cuprates [pO|. 





FIG. 5. Uniform susceptibility xo (in units of 1/t) vs. T 
for various Ch (increasing in steps of 0.05), as obtained by 
calculations on the system with = 20 sites for Ch > 0, and 
N = 26 sites for = 0, respectively. 

To show that such physics is indeed incorporated 
within the t — J model, we present our results for the 
uniform susceptibility Xo(2^) at various hole dopings Ch- 
From Fig. 5 we establish that Xo(2^) is most monotonous 
within the intermediate-doping regime. On the other 
hand, for < 0.15 the susceptibility starts to reveal 
a shoulder at T, which could be interpreted as an in- 
crease of the pseudogap scale, but as well as a sign of the 
dominating role of the AFM fluctuations. 
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0.0 0.1 0.2 0.3 

FIG. 4. Entropy density s vs. for three T and two sys- 
tem sizes A''. 

VII. CONCLUSIONS 

Our results for the t—J model seem to indicate that the 
anomalous, but universal, dynamics in the intermediate- 
doping regime of correlated systems is a consequence of 
an extreme degeneracy and a collapse of the low-lying 
quantum states (introduced by doping the magnetic in- 
sulator). This leads to a diffusive- like charge and spin 
response, where T represents the only relevant energy 
(frequency) scale. On the other hand, our results do not 
exclude the possible onset of coherence (related to 'pseu- 
dogaps' and superconductivity) at lower T < T* ~ 0.1 t. 

It is evident that the 'optimal' doping regime is charac- 
terized by the competing influence of two quantum pro- 
cesses: AFM fluctuations and hole hopping. Outside this 
regime, either one of the two processes should dominate. 
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